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Abstract 

A simple, yet reasonably accurate, analytical technique is proposed for multi-factor structural 
credit portfolio models. The accuracy of the technique is demonstrated by benchmarking against 
Monte Carlo simulations. The approach presented here may be of high interest to practitioners 
looking for transparent, intuitive, easy to implement and high performance credit portfolio model. 

1 Introduction 



Structural multi-factor economic capital (EC) models derived from the CreditMetrics framework (Gupton 
ct al., 1997) have become the most widely adopted tools for risk quantification in credit portfolios. An 
outcome of these models, a portfolio EC and its allocation down to individual facilities, is used by financial 
institutions for any or all of the following: internal capital adequacy assessment, external reporting, risk- 
based pricing, performance management, acquisition/divestiture analyses, stress-testing and scenario 
analysis, etc. While in most cases Monte Carlo simulations are used due to limited analytical tractability 
of the multi-factor models, the recently reported advanced analytical techniques (Voropaev, 2011) may 
fi be viewed as an alternative. 

^^-H Unfortunately, neither the industry standard simulation-based approach nor the existing analytical 

^^ techniques can fully address the needs of the financial institutions. In particular, the risk-based real-time 

1—^ pricing remains the ultimate challenge: none of the existing models is capable of providing sufficiently 

accurate, stable and time-efficient input. Yet another practical aspect which has not received enough 
(^ attention in the literature is the sometimes overly complex structure of the models as perceived by end 

tH- users. Very often the complexity of the models makes them hard to be understood and, hence, affects 

\^ their acceptance within an organization. 

^~~l The approach presented here aims to overcome the above mentioned difficulties and is in its spirit 

^^ similar to the one reported by Ccspcdes ct al. (2006). The content, however, is quite different since the 

C*^ presented model has more solid theoretical background, is easier to implement and use and is capable of 

^^ covering fully featured multi-factor setup. The model described here was developed with KISS principle 

^^ ^ in mind. While based on the previous author's research on the analytical tractability of multi-factor 

. . models (Voropaev, 2011), the proposed model has very simple and intuitive structure. Despite the simple 

. ^ structure, the model produces meaningful and reasonably accurate results and can be used by financial 

k> institutions for any of the purposes described above. In particular, the problem of capital allocation has 

^ a simple and time-efficient solution allowing real-time risk-based pricing. From conceptual point of view, 

^ one of the most attractive features of the model is its ability to quantify risk concentrations on both 

sector and obligor levels in a similar fashion. 

This article is organized as follows. A short description of structural multi-factor model and the 
necessary theoretical background are given in Section 2. Mathematical foundations of the proposed 
model are presented in Section 3 and are substantiated by benchmarking with Monte Carlo simulations 
in Section 4. Section 5 contains some concluding remarks and summarizes the presentation. 
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^Keep It Simple and Straightforward, http://cn.wikipedia.org/wiki/KISS_principle 



2 Background 

Let us consider a portfolio of credit risky facilities with loss functions {h}{ei) at horizon (one year) being 
a function of random variables (normalized asset returns) {si}- Dependencies within the portfolio are 
modeled by means of a set of common factors {rik}'- 



e^^P^Y.mkVk + ^l-pU^P (2-1) 

k 

The random variables {{rjk}, {£,i}} are independent and standard normally distributed ^. The instrument 
specific \pi\ < 1 and {f3ik} define the systematic sensitivities of the instruments. The latter are subject 
to normalization condition '^f. /Sfj, = 1. The idiosyncratic risk components are represented by {£,i}. 

The economic capital of the portfolio is defined as na a - quantile (usually set to 99.9% or higher) of 
the portfolio loss distribution L = J2i h relative to the expected loss of the portfolio: 

EC = qa[L]-E[L] (2.2) 

The above quantifies the overall portfolio risk which can be consistently distributed between the under- 
lying facilities using the Euler principle as (Tasche, 2008, see e.g.): 

EC^^ec, ec, = «;,^EC (2.3) 

i 

where Wi is a weight of the ith asset in the portfolio. To simplify the notations, these weights will not be 
written explicitly in what follows. 

No closed form analytical solution exists for either portfolio EC or its allocation {eci} in general case. 
However, in a single-factor case, i.e. one common factor rjif and pi — 1 for any i, the portfolio loss 
distribution Lif quantile can be trivially found to be'^ qa[Lif] = Lif{r]if = N~^{a)), where N~^{) is 
an inverse cumulative standard distribution function. This provides motivation to look for a solution of 
(2.2) in a form of a sum of a single factor approximation and some corrections (for detailed explanation 
of this approach see e.g. Voropaev (2011) and references therein): 

qa[L] = qa[Lif] + 6qc,[SL„,f], £i/ = E[L\riif], E[L„,f] = (2.4) 

Here Li/ is conditional on the single factor value loss distribution of the portfolio. The single factor 
is constructed as a linear combination of the systematic factors 771/ = ^^ ctkrjk with the normalization 
condition '^^'^k ~ 1- Obviously, the choice of ryi/ significantly affects the quality of the approximation 
(2.4) and will be given a particular attention in what follows. 

The conditional expectation series expansion technique (Voropaev, 2011) can be applied to facilitate 
calculations of Lif{'qif). Translated to the notations introduced here this technique allows writing the 
conditional portfolio loss distribution Lif as a sum of the conditional expectations of the constituents 
U — 'Ei[li\'qif] which can be expressed as: 

Umt) = E ^^^^,^"^He„(,yi,), 4") = / Z,(e)He„(6)^de (2.5) 



n=0 

where the inner vector product /3id stands for ^,^ anl^in and He„() are Hermite polynomials. The series 
converge very well provided values of PiPid are not too close to 1 (which is the case in practice) and allow 
for very fast (re) calculations of the conditional expectations once the constants l^' have been computed. 
The above technique is particularly useful when considering arbitrary loss functions {U}. 

^In practice the common factors are not independent are correspond to industry and geographic sectors. However, their 
correlation matrix can always be diagonalized. The latter is used here to simplify the notations. 
^Here it is assumed that Lif{r]if) is an invertible function. 



3 KISS model 

3.1 Systematic and idiosyncratic risk: happy marriage 

As long as credit portfolio modeling is concerned, it became a common practice to distinguish the un- 
systematic {r]k} and the idiosyncratic {^i} risk components. The usual assumption is that the former 
drive the portfolio risk dynamics while the latter only give minor contributions. However, the two sets 
of random variables do not differ from mathematical point of view. In fact, one cannot draw a clear 
line between the systematic and idiosyncratic components using practical considerations either. Indeed, 
imagine that the portfolio contains a single relatively big exposure or a set of exposures corresponding 
to the same borrower and, hence, sharing the same idiosyncratic random variable ^i. Depending on the 
size of this cxposure(s), the risk brought to the portfolio by ^i may be higher than the one originating 
from some or even all of the systematic factors {ilk}- Big enough exposure will eventually dominate the 
portfolio dynamics even if its sensitivity pi to the systematic factors is zero. Introducing credit contagion 
effects by assigning more than one overlapping idiosyncratic factors to a group of dependent borrowers 
makes it even harder to make a distinction between the systematic and idiosyncratic factors. 

Treating the systematic and the idiosyncratic risk equally not only simplifies the model structure, 
but also allows straightforward incorporation of the borrower concentration effects into the portfolio risk 
metrics. The notations used so far can be generalized as follows. For M common factors {rjk} and the 
portfolio consisting of N borrowers let us introduce 

pl = (pi/3ji,PiA2,---,A/3jAf,0, ...,yi-pf,...,0), |lPi|l=l (3-1) 

a = (ai,. . . ,aM,o:M+i,- ■ ■ ,aM+N), l|a|| = 1 (3-2) 

M+N 

Ti = pj . a = ^ PikOik (3.3) 

The single factor approximation Lif of the portfolio loss can be written as 



oo 

r. 



Li^.=^Z„ I,(r„,yi/) = ^^/,|")He„(7yi/) (3.4) 



n 

n=0 



Unification of the systematic and idiosyncratic risk factors results in idiosyncratic factors being in- 
corporated in a which, as will be shown, defines the portfolio risk dynamics. Thus, the idiosyncratic risk 
is accounted for in the same fashion as the systematic one. 

3.2 Best single factor approximation 

As was mentioned before, analytical tractability of the single factor case is the starting point for approach- 
ing the more general multi-factor setup. Starting with a single factor approximation and calculating 
multi-factor (including idiosyncratic) adjustments as in (2.4), one can in principle calculate the portfolio 
economic capital. This approach, however, suffers from some difficulties. First, the choice of rjif is not 
obvious, yet a very important first step. Next, calculations of the multi-factor corrections may be quite 
laborious and hardware demanding. Finally, as will be demonstrated later, the multi-factor corrections 
are not guaranteed to be convergent. 

Instead of trying to overcome the difficulties associated with the multi-factor adjustments calculations, 
let us put all the efforts into constructing the single- factor approximation. Some factor 1/ should exist 
which maximizes the relative contribution of the single-factor approximation in (2.4) and, thus, dimin- 
ishing the relative importance of the multi-factor corrections. Assuming that the multi-factor corrections 
give positive contribution"' to the a-quantile of the portfolio loss distribution L, the optimization problem 
reduces to maximization of the contribution from the single- factor approximation ii/: 

qa[L] w iRaxqa[Lif] (3.5) 



^This is not a solid assumption from mathematical point of view; however, it is true for any portfolio one may encounter 
in practice. In the worst case scenario, i.e. negative multi-factor corrections being neglected and the maximum single-factor 
contribution being used, one is only risking being somewhat conservative. 



The validity of this crucial assumption will be substantiated later in Section 4 by benchmarking with 
Monte Carlo simulations. 

From now on let us define the economic capital EC of the credit portfolio as an a-quantile of the 
optimal single factor distribution Lif. Using the notations introduced in this section the economic 
capital can be written as 

EC = Vl,(r,) -Ve[Z,], r, = a-pi (3.6) 

The optimal single factor is defined by a which maximizes the above expression 

V5^I,(r,;)=0 (3.7) 



and has the following solution 



IIpII ^-^ on 



This equation, however, contains a on both sides (r^ on the right contains a) and does not allow a 
straightforward analytical solution. The problem (3.5) can still be solved numerically by applying, for 
example, the method of steepest descent. Based on (3.8), the following starting point can be suggested 



-^ Po - V^ dliin) ^ 

bo "^ dr. 



(3.9) 



The calculations can be significantly facilitated by the series expansion (3.4). In practice, only few 
iterations are needed to have an accurate solution to the optimization problem (3. -5). The calculations 
are not hardware demanding and very fast. 

The optimal single factor defined by (3.8) leads to another simplification for the portfolio capital 
allocation problem (2.3). The individual capital contributions 



OWi OWi ^—' 



ec,; = z«,,^^EC = w~Y^ (Un) - E[k]) (3.10) 



can be written as 



ec, = h - E[h] + y -^pj .w,j—^^ hidpi) - E[h] (3.11) 

'^ dri ' dw, \\p\\ 



where the third term can be shown to be zero: 



dlj _, ^^ d P _ -, „„ d P _ ^ ( Pi pip- Pi) 



^ dr-i dw^\\p\\ ^w^\\p\\ \\\p\\ \\p\" 



In other words, the choice of a according to (3.8) leads to particularly simple expressions for capital 
contributions (3.11). The overall portfolio EC is a sum of conditional expectations of the excess losses 
hivif = ^~^i'^)) ~ ^[h] which happen to coincide with the individual capital contributions. 

Let us emphasize that once the optimal single factor (3.8) has been computed, the set of parameters 
a is sufficient to perform the capital allocation calculations (3.11). This allows real time calculations 
necessary for risk-based pricing. In fact, the proposed method for the economic capital calculations and 
its allocation is so efficient that the calculations can be performed on huge portfolios containing millions 
of facilities in minutes using entry level desktop computer. The accuracy of the method is demonstrated 
in the next section through benchmarking against Monte Carlo simulations. 



4 Benchmarking 

To highlight the points made and substantiate the assumptions used in the previous section, let us compare 
the performance (accuracy) of the proposed analytical approximation with Monte Carlo simulations. 
To demonstrate the advantages of the proposed technique, the comparison analysis will also cover the 
previously reported analytical technique (Voropacv, 2011) which, in contrast with the one presented here, 
aims at precise calculations of the multi-factor corrections in (2.4). 
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Figure 1: Portfolio EC: Monte Carlo vs. analytical estimates. 



4.1 Artificial portfolio, high concentrations 

Let us start the analysis by considering a simple portfolio of 1000 bullet loans maturing at the horizon. 
Each loan has sensitivity p| — 0.2 to the single systematic factor rj and a unique idiosyncratic component 
^i. The loss at horizon functions are 



kiei) 



lo if e> -iV-i(PD) 
if e< -iV-i(PD) 



(4.1) 



where PD; are probabilities of default and are equal for all the loans. Two sets of experiments were 
conducted with PD=1% and PD=0.1%. The loss severities Iq are initially set equal to all loans in the 



portfolio. The concentration effects are studied by gradually increasing the loss severity of one of the 
loans and examining the impact on the portfolio EC. The confidence level for EC is set to 99.9% except 
for the last experiment where the dependency of the EC on the confidence level was investigated for fixed 
f 0% concentration. Both EC and concentration are measured as a fraction of the total portfolio exposure 
(i.e. the sum of all the loss severities). 

The results of the Monte Carlo simulations are compared with the numbers produced using the KISS 
model (Section 3) and with the output of the approach based on the calculations of the multi-factor 
corrections in (2.4). In the latter case the single systematic factor rj was used for the single-factor 
approximation. 

Based on the results of the comparison presented in Fig. f , one can conclude the following. Multi- 
factor corrections may lead to very accurate results in case of moderate concentrations (PD=1%, concen- 
tration < 10%). In case of high concentrations, however, the results suggest that the higher order (3rd 
order in this case) multi- factor corrections are divergent, while limiting the multi- factor contributions to 
the 2nd order corrections leads to significant underestimation of EC. 

The KISS model, on the other hand, despite being not extremely accurate, produces robust and 
reliable EC estimates for a wide range of parameters. This is especially obvious when studying the 
dependency of the portfolio EC on the confidence level. 

4.2 Realistic portfolio, moderate concentrations 

Here the analysis of the previous section is complemented by the one conducted on more realistic portfolio. 
The portfolio consisted of 2,000 loans to distinct customers randomly selected from a loan portfolio of 
a large European bank'. The set of common systematic factors covering 45 geographic regions and 61 
regions, as well as the valuation function at horizon li{e) used in the experiment were similar to those of 
the PortfolioManager (Kcalhofcr, 2001) model. 

Both the portfolio EC and its allocation {eCi} were estimated using unbiased Monte Carlo simulations. 
The confidence level was set to 99.9% and the EC contributions were estimated as average realized values 
in the interval 99.85%-99.95%. A total of 10^° (ten billion) scenarios were used. 

The simulation-based estimates were compared to both the KISS model outcome as well as with the 
output of the approach based on the calculations of the multi-factor corrections in (2.4). In the latter case 
the single-factor 1/ used as a starting point was obtained using (3.9) restricted to systematic components. 
Second (m/2) and third order (m/3) corrections were calculated. 

The results on the portfolio level are summarized in Table 1 and the facility-level results are presented 
as scatter plots in Fig. 2. 
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Table 1: Relative differences between analytical and simulation- 
based estimates of the portfolio EC. 

On the portfolio level the accuracy of the proposed approximation is more than sufficient for any 
practical purposes. The performance may seem to be less impressive on the facility level. As expected, 
the most significant mismatch is observed for the facilities with high concentrations of capital. These 
discrepancies, however, do not jeopardize the practical validity of the KISS model. Indeed, the uncertainty 
in the input parameters observed in practice (i.e. default rates, recoveries, correlation parameters, etc) 
as well as dependencies on the modeling assumptions diminish to large extent the approximation errors 
of the proposed technique. This, combined with the robustness demonstrated in the previous section, 
makes the KISS approximation a valid alternative for quantification of risks in credit portfolios. 

^The same portfolio was used by Voropacv (2011). 
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Figure 2: Relative differences between Monte Carlo and analytical estimates. 
Ten biggest consumers of capital accounting for 19% of EC are marked. 

5 Summary 

Despite its simplicity, the analytical approximation presented here is capable of quantifying credit port- 
folio risks in a general multi-factor setup. The VaR risk measure used here can easily be replaced with 
the Expected Shortfall. The arbitrary loss functions {li{f-i)} used allow for covering not only default-only 
regime, but also MtM valuation or even the dependency of in-default loss severities on the systematic 
factors. The default-only case has a particularly simple solution mimicing the well-known IRB capital 
rules: 



eci = EaDi • LGDi 






(5.1) 



The less than perfect accuracy of the approximation is not crucial for day-to-day practical needs of 
credit portfolio managers. The advantages of the proposed technique are significant, the model allows 
very fast and straightforward calculations including real-time risk-based pricing. The simple, robust and 
transparent structure can facilitate user acceptance and integration on all levels of financial institutions. 
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